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Abstract 

The 2 £+ — 2 n separation in KO is investigated using large basis sets and high 
levels of correlation treatment. Relativistic effects are included at the Dirac-Focx 
level and reduce the separation only slightly. The basis set superposition error is 
considered in detail. On the basis of these calculations, our best estimate places 
the 2 n 3/2 state about 200 cm -1 above the ground 2 S - state in agreement with our 
previous estimate. 

I. Introduction 

In 1986 we computed the 2 S" r — "II separation in KO to be 238 cm with 
the 2 £ + state being the lower of the two [1]. This separation was computed at 
the singles and doubles configuration interaction (SDCI) level using a large Slater 
basis set. While not definitive, the separation was thought to be accurate because 
of the small difference between the SDCI and self-consistent-field (SCF) values. 

However, when 15 electrons are correlated, an SDCI treatment can underestimate 
the differential correlation effects. Approaches such as the modified coupled-pair 
functional (MCPF) method [2], or the coupled cluster singles and doubles method 
with a perturbational estimate of the triple excitations [3] [CCSD(T)], account for 
the effect of higher than double excitations and can therefore yield superior results 
to the SDCI approach. 

Experimentally, the ground state has not been determined and there is conflict- 
ing evidence. The failure to observe an electron spin resonance (ESR) spectrum [4] 
is consistent with a 2 n ground state, while the magnetic deflection experiments [5] 
suggest a 2 £ + ground state. Photoelectron detachment experiments [6] on KO - 

! Mailing address: NASA Ames Research Center. Moffett Field, CA 94035-1000 

(i\«ASA-TK-10A0«o) TH tv 2 5 ( + ) - IP N93-1 2°00 

j t P A A T I U f'i IN VI (NASA) 9 ;> 

U n c 1 a s 

G 3 / 1 Z 


C 1 c 1 J 5 6 



might offer a mechanism to directly observe both states, and hence determine the 
ground state and the separation; such experiments [6] have recently been performed 
for NaO~ and preliminary results suggest that they are possible for all of the alkali 
oxides systems. In addition, microwave experiments [7] have recently been per- 
formed on NaO and, with a pure-precession hypothesis, an estimate of the £ — II 

separation was obtained. An analogous experiment should be possible for KO. 

Given the recent experimental activity on the alkali oxides, it seems timely to 
reinvestigate the 2 £" J ' — 2 II separation in KO at higher levels of theory. We consider 
the effect of one-particle basis set and basis set superposition error (BSSE) in more 
detail than in our previous work. In addition, the separation between the £ and 
2 II states is small enough that relativistic effects could affect the ordering of the 
states. These were not included in the previous work and are addressed here. 

II. Methods 

Two basis sets are used in this work. The small oxygen basis set is the 
(135 9p 6d 4/) /[5s 5p (2 -f 1 )d If] atomic natural orbital [8] (ANO) set [9] de- 
termined from the average natural orbitals of 0 and O . (The notation ti — m in- 
dicates n ANO orbitals are employed and m primitive functions are uncontracted). 
In the large oxygen basis set only two s and one p ANOs are used. The outer 
five s and five p primitives and the d and f polarization functions are uncon- 
tracted. This is augmented with two uncontracted g functions and a diffuse s 
function with an orbital exponent of 0.0(6666. Thus the final basis set is of the 
form (14s 9 p 6 d 4/ 2p)/[(2 + 6)s (1 4- 5)p 6d 4/ 2g\. 

The small K basis is derived from the (24s 15p) basis [10] augmented with an 
_g function optimized for K and three p functions optimized for the “P state of K. 
The s and p basis sets are contracted to (3+6) and (2+7), respectively, where the 
three s and two p contraction coefficients are taken from the SCF wave function. 
The basis is further augmented with six even-tempered (a n =2.5 n a 0 for n=0,fc) d 
functions with o:o=0.0455 and four even-tempered f functions with uo — 0.12949. 
The small K basis set is of the form (25s 18p 6 d 4/)/[(3 + 6)s (2 + 7)p 6 d 4/]. In 
the large K basis set, the s and p functions are contracted (2+12) and (2+9) and 
two g functions (o;=0.768 and 0.307) are added, thus yielding a basis set of the form 
(25s 18p 6d 4/ 2g)/[{2 + 12)s (2 + 9)p 6d 4/ 2 g). Only the pure spherical harmonic 
components of the basis functions axe used. 
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The orbitals are optimized using the SCF approach. We have treated the II 
state using symmetry and equivalence restrictions (i.e., a occupation ) or 

by using a symmetry broken solution (i.e., a 7r"7r* occupation). Electron correlation 
is added using the MCPF or CCSD(T) approaches. Fifteen electrons are correlated, 
corresponding to the 0 2s and 2 p and the K 3-s, 3p, and 4s electrons. 

Relativistic effects are estimated as the difference between the SCF and Dirac- 
Fock (DF) results. Because correlation was not included, it was possible to reduce 
the basis set size without affecting the accuracy of the results. The full s and p 
primitive basis sets for K and 0 were contracted to [(3 -f* 6)s (2 — *)pj and [(2 -r 
3)s (1 + 4)pj using the orbitals from the SCF wave functions for K~ and 0~. For 
polarization the 6 d set and the outer three of the 4 f functions on R. were used, and 
on 0 the outer 3 primitive d functions and the / function with the largest coefficient 
from the ANO polarization set were used. The SCF 2 E + — 'II separation in this 
basis set was identical to that found using the small ANO set. 

The SCF/MCPF calculations were performed using the SEWARD [11]- 
SWEDEN [12] program system. The CCSD(T) calculations were performed using 
the program developed by Scuseria [13]. The atomic DF calculations were performed 
using GRASP [14] while the molecular calculations were performed using the code 
developed by Dyall [15]. All calculation were performed on the NASA Ames Central 
Computer Facility CRAY Y-MP or Computational Chemistry Branch IBM RISC 
System/6000 computers. 

III. Results and Discussion 

The separation at the SCF level is only slightly larger than in our previous 
work [1], with the big basis set result essentially reaching our previous estimate of 
250 cm -1 for the numerical Hartree-Fock separation. The SDCI separation in the 
small basis set is only 6 cm -1 smaller than that reported in our previous study. 
With a symmetry broken SCF solution, the 2 II state is predicted to be below the 
2 E+ state. When correlation is added at the SDCI level, the 2 E + state is predicted 
to be lower. However, the SDCI based on the symmetry broken SCF orbitals yields 
a significantly smaller separation. This suggests that the SDCI approach is unable 
to overcome the bias in the SCF wave function. At the MCPF level, the effect of 
sy mm etry breaking is smaller, as expected, since the MCPF approach includes the 
effect of higher excitations. At this level, the symmetry broken solution actually 
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leads to a larger separation than using the SCF with symmetry and equivalence 
restrictions. Both MCPF separations are about 100 cm -1 smaller than the SDCI 
result with symmetry and equivalence restrictions. The CCSD separation is the 
same as the MCPF. The effect of triples is to further increase the separation by 
47 cm -1 , which is almost an order of magnitude smaller than the difference be- 
tween the SCF and CCSD (368 cm -1 ). This strongly suggests that even higher 
levels of theory will not significantly increase the correlation contribution to the 
separation. Expanding the basis set results in a small increase in both SCF and 
MCPF separations. 

In Table II we report the BSSE at the MCPF level. We compute the BSSE for 
K + with the 0 basis set, K + (0), and O - with the K basis set, because both states 
of the molecule at r e are best described as K + 0~. K - is a closed shell; therefore the 
calculation of the BSSE is straightforward, and we consider K -r (0) at the two rele- 
vant internuclear separations. The 2 p 5 occupation of 0 is oriented 2p<r"2p7r 3 and 
2 p< 7 1 2 p 7 r 4 for the 2 II and 2 E + states, respectively. As shown in previous work [16], 
the BSSE can differ significantly for the two occupations. Therefore we use the 
BSSE for the 2 pcr 2 2 p 7 r 3 occupation to correct the 2 II state and the BSSE for the 
2p<x 1 2p7r 4 occupation to correct the state. 

For K + the BSSE is small for both basis sets. It is slightly larger for the 2 E+ 
state than the 2 II state because of the shorter bond length. In contrast, the O - 
BSSE is much larger for the smaller basis set. In spite of the longer bond length, 
the 0 - (K) BSSE is larger for the 2 1I state, because a pair of oxygen 2p<r electrons 
points at the ghost basis set compared with only one electron in the "S • state. 
For both basis sets the total BSSE is very similar for the "E ■ and "II states and 
therefore does not have an important differential effect on the separation. Given 
the small total BSSE for the big basis set, this conclusion seems definitive in spite 
of any uncertainties in how to best compute the BSSE. 

One effect not considered to this point is relativity. The spin-orbit splitting in 
the O - ion was calculated to be 193 cm -1 at the Dirac-Coulomb level. Inclusion 
of the Breit interaction decreased the splitting to 170 cm -1 . This value is expected 
to be accurate as an analogous treatment of O agrees with experiment to within 
2 cm -1 . Because the bonding in KO is best described as K + 0“, this splitting in O - 
can be used in a simple first-order model [17] to estimate the splitting in KO. On 
this basis, the 2 E^ 0 state is unshifted from its nonrelativistic position, the Hi / 2 
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state is shifted up by 57 cm -1 and the 2 n 3 / 2 state is shifted down by 57 cm 1 . 

Mixing of K 4s character from covalent contributions to the wave functions 
and spin-orbit mixing of the 2 Hi yo and 2 S^, 2 states could invalidate the first-order 
model. Therefore we study the "II 3/2 and states at the Dirac-Fock level. 

On the basis of the nonrelativistic calculations we use a bond length of 4.4 do for 
the 2 n 3/2 state and 4.13 a 0 for the 2 E+ /2 state. At this level, the 2 n 3 / 2 state 
is lowered by 36 cm -1 and the 2 S+ /2 state is lowered by 24 cm -1 , relative to the 
nonrelativistic states. Thus the differential effect of relativity on the state separation 
is only 12 cm -1 , or about 1/5 of the estimate derived from 0“. Clearly, however, 
the differential effect of relativity is sufficiently small that it will not reverse the 
order of the states computed at the nonrelativistic level. 

Our best separation is obtained by correcting the CCSD(T) separation for: 1) 
MCPF basis set improvement (+36 cm -1 ), 2) BSSE (+8 cm ■*) and 3) relativistic 
effects (-12 cm -1 ). This yields a separation of 214 cm -1 , with the 2 E + state 
being lower. This is very similar to our previous result [lj. While it is difficult 
to be definitive for a separation this small, the calculations strongly suggest that 
the ground state is ‘S’’’. Because improving the basis set increases the separation, 
as does including higher levels of correlation treatment, and relativistic effects are 
very small, it is difficult to imagine what effect would reverse the order of these two 
states. 

IV. Conclusions 

The separation at our highest level of theory, the CCSD(T) approach, is 
182 cm -1 . Including the correction for basis set incompleteness from the MCPF 
level increases the separation to 218 cm -1 , which is very similar to our previous 
result [1]. The differential effect of BSSE is 8 cm -1 , increasing the separation. Rel- 
ativistic effects reduce the separation by only 12 cm 1 . Thus we conclude that the 
ground state is 2 E + . 
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Table I. Spectroscopic constants for the and 2 1I states of KO. 
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Table II. Summary of the BSSE at the MCPF level. 


BSSE(cm _1 ) 

Small basis set 

2 E + 

K+(0) at r=4.13 32.1 

0~(K) at r=4.13 348.4 

Total 380.5 

2 n 

K+(0) at r=4.40 20.9 

0~(K) at r=4.40 364.5 

Total 385.4 

Big basis set 

2 £+ 

K + (0) at r=4.13 32.5 

0"(K) at r=4.13 88.5 

Total 121-0 

2 n 

K + (0) at r=4.40 25.8 

0~(K) at r=4.40 103.4 

Total 129.2 
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